Non-monotonic magnetoresistance of two-dimensional electron systems in the ballistic 

regime 
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We report experimental observations of a novel magnetoresistance (MR) behavior of two- 
dimensional electron systems in perpendicular magnetic field in the ballistic regime, for ksTr/h > 1. 
The MR grows with field and exhibits a maximum at fields B > 1/fi, where fi is the electron mobil- 
ity. As temperature increases the magnitude of the maximum grows and its position moves to higher 
fields. This effect is universal: it is observed in various Si- and GaAs- based two-dimensional elec- 
tron systems. We compared our data with recent theory based on the Kohn anomaly modification 
in magnetic field, and found qualitative similarities and discrepancies. 

PACS numbers: 73.63.Hs, 73.40.Qv, 73.40.Kp, 73.23 
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Two-dimensional (2D) degenerate electronic systems of 
high purity (kpl 3> 1) with isotropic energy spectrum are 
rather simple objects, which seem to be well understood. 
Within the classical kinetic theory^ the resistivity of such 
a system should not depend on perpendicular magnetic 
field for lo c t < 1 (where us c — eB/m* is the cyclotron 
frequency, and r - the transport time). However, a no- 
ticeable magnetoresistance (MR) is often seen in exper- 
iments with 2D systems; such MR is usually attributed 
to quantum corrections which are beyond the classical 
consideration. There are two types of quantum correc- 
tions to conductivity: (i) weak localization (WL), and 
(ii) electron-electron (e-e) interaction correction (for a 
review, see Reb 2 -). In the diffusive regime (kBTr/H <C 1, 
t/t v <C 1), both corrections give rise to magnetoresis- 
tance with an amplitude proportional to \n(T]2£ whereas 
in the ballistic regime {ksTr/h > 1, t/t v > 1) the 
magnetoresistance should disappear—. These theoreti- 
cal predictions for the MR have been verified in diffu- 
sive and diffusive-to-ballistic crossover regimes in recent 
experiments^^ with 2D systems. A conventional belief 
(that the quantum corrections to MR have to disappear 
at high temperatures) has made the MR in purely bal- 
listic regime out of the scope of experimental interests. 
This theory prediction for the ballistic regime, however, 
was not verified thoroughly. In order to shed light on this 
issue, we measured MR in the ballistic regime with var- 
ious simple isotropic 2D electron systems. Contrary to 
the common belief we have found that the MR in perpen- 
dicular fields does not vanish at ksTr/h > 1; instead, it 
manifests a novel type of behavior: the MR depends non- 
monotonically on field and exhibits a maximum, whose 
position scales with temperature for all samples. 

In this paper, we report observation and systematic 
studies of the MR in the domain kpl ^S> 1, /csTr/fi > 1, 
where the MR should be missing. Experimentally, how- 
ever, different Si-MOS structures, GaAs/AlGaAs het- 
erostructures and GaAs-based quantum wells were found 
to show a nonmonotonic MR. Our results provide an ev- 
idence for a universal origin of the effect. We compared 



our data with a recent theory! of e-e interaction correc- 
tion that employs modification of the Kohn anomaly by 
magnetic field and did find some qualitative similarities. 

We used two Si-MOS samples (Si4, Sil3 with peak 
mobilities 1-2 m 2 /Vs) and GaAs-AlGaAs heterostruc- 
ture 28, GaAs24 (mobility 24 m 2 /Vs) 8 , and gated 
quantum well structures AlGaAs-GaAs-AlGaAs (1520) 
and GaAs-InGaAs-GaAs (3513)4 All samples were 
patterned as Hall bars. Density of electrons in gated 
samples was varied in situ. The relevant parameters of 
the samples, densities n (in units of 10 12 cm~ 2 ), and mo- 
bilities /i (m 2 /Vs), are summarized in the following table: 
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Samples were inserted into a cryostat with a supercon- 
ducting magnet; the field direction was always perpen- 
dicular to the 2D sample plane. Temperature was varied 
in the range 1.3-60 K. Both components of the resistivity 
tensor were measured simultaneously using the standard 
four-terminal technique with either SR-830 lock-in ampli- 
fier (samples Si4, Sil3, 28, 24), or using rectangular cur- 
rent modulation (samples 1520,3513). Both, harmonic 
and rectangular modulation was made at frequencies 12- 
33 Hz. Current was chosen an order of 1 fiA, to ensure 
the absence of electron overheating. 

In order to exclude an admixture of the off-diagonal 
component of the resistivity, we swept magnetic field 
from — B to B, and then symmetrized our data. Such 
a symmetrization is necessary for reliable measurements 
of corrections to the resistivity whose relative variations 
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FIG. 1: (Color online) Magnetoresistance for sample Si4 
at different temperatures. Electron density n = 1.72 ■ 
10 12 cm -2 . Up-arrows mark positions of the p xx maxima. 
Ap designates the magnitude of the MR. Dashed curve 
on the panel e shows fitting according to Eq. ((2)1 with 
A 2 = 0.2. h/k B r « 8K. 
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FIG. 2: (Color online) Magnetoresistance for sample 1520 
at different temperatures. Electron density n = 1.4 ■ 
10 12 cm" 2 . Temperature values are indicated in the fig- 
ure, h/ksr = 13. 5K. 



might be less than 1%. 

Electron density values quoted in the paper were de- 
termined from the slope of the Hall resistance versus B 
as well as from the period of Shubnikov - de Haas oscilla- 
tions at low temperatures. Both results agreed with each 
other within 2%. The highest temperature in our experi- 
ments was chosen not to exceed 60 K in order the carrier 
density to remain constant and to avoid a bypassing bulk 
conductivity. 

Examples of our MR curves, obtained at different tem- 



peratures for samples Si-4 and 1520 at fixed electron den- 
sities are shown in Figs. [T]and|2J respectively. As mag- 
netic field is increased from zero, at first, p xx sharply falls 
due to weak localization suppression. Then p xx starts 
growing and reaches a maximum at iJ max field (indicated 
by the arrows in Fig. [1} . After passing the maximum p xx 
decreases; in higher fields, \B\ > 1.5B max , MR can be- 
come either positive or negative depending on the sample, 
temperature, electron density, etc. At the lowest temper- 
atures, Shubnikov - de Haas oscillations are seen in high 
fields, on top of the smooth MR. 

This unexpected nonmonotonic magnetoresistance is 
the main subject of the current paper. We stress that 
this effect (i.e., nonmonotonic MR) is universal. The 
point is that in different samples and at various electron 
densities it has similar features: (i) MR is small (its typ- 
ical magnitude is less than 1%), (ii) the nonmonotonic 
MR is observed only for T > 1.3h/k B T 9 , (iii) the MR 
maximum grows in magnitude and moves to higher mag- 
netic fields as temperature increases (the position of the 
maximum exceeds lu c t > 1 and is roughly proportional 
toT). 

Comparing the data from Figs. {1} and Fig. (JSJ) for Si- 
MOSFET and GaAs QW-samples with similar mobilities 
and densities, we see that the MR takes a maximum at 
similar temperatures and magnetic fields, and at similar 
u! c t values. This result indicates that the MR has an 
orbital rather than spin origin because the Zeeman en- 
ergies g* fis differ by a factor of 5 for these two different 
material systems. Also, this effect has nothing to do with 
WL and e-e-interaction diffusive corrections^ because it 
survives at such high temperatures as huTr/h w 20 for 
samples 28 and 24 at T = 20K. 

Searching for possible semiclassical mechanisms, we 
have to note that most of the theoretical models for 
the case of short-range scatterers^ predict a negative, 
monotonic and temperature independent magnetoresis- 
tance, due to the memory effects^. A positive, though T- 
independent, magnetoresistance was predicted in RefJ^, 
due to non-markovian scattering. The latter type of MR 
was experimentally observed in very clean samples and 



for classically large magnetic field: 
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~t > 1. There- 



fore, we conclude that the aforecited semiclassical mech- 
anisms can't explain the nonmonotonic MR observed in 
our experiments. 

Recently, Sedrakyan and Raikh- suggested a new MR 
mechanism, which causes a maximum of resistivity in 
not-too-strong magnetic fields oj c t ~ 1 in the ballistic 
regime (huTr/h > 1). This new mechanism seems to 
give the best starting point for comparison with our mea- 
surements. The MR in Rett originates from the e-e inter- 
action correction to conductivity. According to RefJ^, e- 
e interaction corrections to conductivity arise from scat- 
tering of electrons on Fridel's oscillations of electron den- 
sity around impurities. Fridel's oscillations are a mani- 
festation of the Kohn 2fcp anomaly in screening. Mag- 
netic field applied perpendicular to the 2D plane modifies 
the electron spectrum and the Kohn anomaly; hence, the 
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field affects screening and electron scattering. In Ref. 7 
this point was taken into account and shown to give rise 
to the second-order correction in the ballistic regime (see 
Eq. (5) from RefX): 



5<7 rr 



4A 2 



irkfiT 



E* 



3/2 



Fo 



uj r E 



1/2 



(1) 



1 



3F V(1 



I 27r 3 / 2 (fc B T) 3 / 2 J 
Fq) is the interaction 



where A = 
parameter—. 

Several predictions can be made based on this equa- 
tion: (1) the correction to resistivity in small fields is 

always positive, (2) (5a xx /a xx ) ■ (Ep/T) 3 ^ 2 is a universal 
function of lu c E^/ 2 /T 3 / 2 for a given interaction strength 
A, and (3) MR has a maximum at uj c t rj 1/ \/3. 
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FIG. 3: (a) oj™ ax r value versus dimensionless temperature 
ksTr/h for all samples. Electron densities (in units of 
10 12 cm" 2 ) are - n = 0.6 (Sil3);4 - n = 0.7 (Sil3);x 
n = 1.3 (Si4); H - n = 1.7 (Si4); EB - 
n = 3.4 (Si4); - n = 1.4 (1520); ® 
ra = 1.6 (1520); - n = 0.8 (1520);* 
1 (1520); V - n = 0.35 (28); A - n = 0.4 (24). 
Dashed line corresponds to ftu;™ ax = 0.7fcsT\ Dotted line 
is ?ia;™ ax = 0.2fcsT. Horizontal thick line is the theoretical 
prediction (see in the text), (b) The same data for GaAs- 
based samples solely scale in coordinates cj" lax T versus 
k B T/E F . 

By comparing these theoretical predictions with our 
data, we find that prediction (1) is always fulfilled af- 
ter subtraction of the weak localization. As for predic- 
tion (2), the p xx (B)-d&t& for different temperatures and 
over the whole range of magnetic fields do not scale as 
the theory predicts. Furthermore, the position of the 
MR maximum in our data is temperature dependent and 
corresponds to uj c t w 1 — 3, contrary to prediction (3). 
Moreover, the magnitude of the MR falls as tempera- 
ture raises in the theory, whereas in our experiment it 



grows with temperature. Evidently, there is no complete 
agreement between the theory! and our experiment. We 
note finally, that the magnitude of the MR maximum Ap 
seems do not scale with any dimensionless combination 
of ksT, H/t, ftw™ ax and Ep. This is also in contrast with 
the theory, where Ap should be cx ft 2 /[r 2 £;^ 5 (fc B T) 1 - 5 ]. 

According to EqQ] the magnitude of the effect is pro- 
portional to the interaction constant A 2 . Therefore, one 
could estimate A 2 value from the experimental data. In 
the theory, the maximum of MR inevitably results from 
(1 — w 2 r 2 ) prefactor in resistivity tensor and should occur 
at uj c t « 1/V3. On the other hand, in the experiment 
the p xx maximum is always observed at lo c t > 1, which 
indicates that this prefactor is weaker than in the the- 
ory. Therefore, for the order-of-magnitude comparison, 
we rewrite Eq. ([T]) for resistivity and omit the [1 — (lu c t) 2 ] 
prefactor: 
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Example of the corresponding fitting with a single vari- 
able parameter A 2 is shown in Fig. Q]e. The fit was per- 
formed in the limited range of magnetic fields 0.15w™ ax < 
uj c < 0.65w" lax , i.e. in the range of the applicability of 
Eq. @ which ignores weak localization and the MR max- 
imum. The A 2 values obtained from the fit appeared to be 
temperature dependent i.e. grew monotonically from 0.1- 
0.4 to 1-3 as temperature was increased from 1.3S/(fcsr) 
to maximal temperature. Surely this temperature depen- 
dence causes the lack of the scaling predicted by Eq. |T]). 
Moreover, A 2 -values obtained from the fitting don't show 
a systematic dependence on carrier density and on ma- 
terial system. 

On the other hand, the A 2 value in our range of den- 
sities may be evaluated from the earlier measurements 
of Fq(u) parameter. The calculated A values are T- 
independent and lie in the interval from 0.2 to 0.5 for 
GaAs-based structures^ and from 1.5 to 5 for Si-based 
structures^. We conclude therefore that the observed 
MR disagrees qualitatively with the theory, though the 
theory predicts the MR of the right order of magnitude. 

In Fig. [3]a, the position of the MR maximum of tem- 
perature. The cj™ ax T value systematically exceeds the 
theoretical expectation 1 / y/3 (horizontal thick line in Fig. 
[3]a) and approximately equals 0.7fcsTr//i for most of the 
data (dashed curve in Fig. [5^,). For samples with the 
highest mobility (24,28), the slope w™ ax r/ '(fc B Tr /K) w 
0.2 whereas for GaAs-based sample with the lowest mo- 
bility the slope exceeds 0.7. In order to take this fact 
into account we have applied another scaling, in coor- 
dinates versus ksT/Ep (see Fig. [3p). It is remarkable, 
that for GaAs-based 2D systems with mobilities and con- 
ductivities ranging by more than an order of magnitude, 
the o;™ ax r data indeed scale reasonably, the result that 
might suggest a clue for understanding the effect. 

The data for Si-based structures are not shown in 
FigEb because they fall out of the T/Ep scaling. In 
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order to understand the origin of the difference in scaling 
for Si- and GaAs- samples, we note that for GaAs-based 
samples in high fields, B > B max , the MR is always neg- 
ative while for Si-based samples it can be either negative 
or positive, depending on particular sample and electron 
density. It means that some other mechanisms affect MR 
in Si-MOSFETs in strong perpendicular fields B > B max 
and shift the MR maximum. It is also worthy of noting 
that in Si the discussed weak MR is observed at such high 
temperatures where metallic temperature dependence of 
the resistivity is strong and nonlinear with respect to T, 
and hence, the first order interaction corrections^ are 
inapplicable. 

We note also, that due to clear reasons the nonmono- 
tonic MR in the ballistic regime T > l.Sh/ksT was not 
oserved in the following cases: (i) Si-MOSFETs in the do- 
main of strong interactions (n < 6 • 10 n cm -2 ) where the 
giant negative MR developes and dominates over other 
weak effects^, (ii) Si-MOSFETs for such high tempera- 
tures where Fermi-gas is non-degenerate (T/Ep > 0.5), 
(iii) GaAs based samples at such high temperatures that 
the carrier density becomes B- and T-dependent. 

Conclusions. In this paper we report experimental 
observation of the novel non-monotonic behavior of the 
magnetoresistance for 2D electron systems in perpendic- 
ular field. This MR is intrinsic to various 2D systems 
(Si-MOSFETs, GaAs and InGaAs quantum wells, and 
GaAs/AlGaAs-heterostructures) and occurs in the ballis- 
tic regime of high temperatures TV > 1. The MR is pos- 



itive in low fields and reaches a maximum at lo c t = 1 — 3; 
the position of the maximum w™ ax scales linearly with 
temperature for all samples. We compare our data with 
recently suggested MR mechanism^ and find some simi- 
larities: (i) the MR is always positive in low field, (ii) the 
MR exhibits a maximum in higher field and (iii) the MR 
magnitude is of the same order of magnitude as predicted. 
However, other features of our experimental data are in 
discrepancy with the theory Refi (i) the MR maximum 
is achieved in fields which are noticeably higher than pre- 
dicted, (ii) the position of the MR maximum linearly de- 
pends on temperature rather than remains constant, (iii) 
the magnitude of the effect increases with temperature 
rather than decreases, as predicted. 

Some clue to understanding the effect may be provided 
by scaling of the MR maximum position w™ ax r versus 
T/Tp, which is empirically observed for various GaAs- 
samples in wide ranges of temperature, density and mo- 
bility. The observation of the nonmonotonic MR shows 
that the magnetotransport theory is still incomplete, at 
least for the ballistic regime, and requires further consid- 
eration. 

Acknowledgements. We thank M. E. Raikh, 
T. A. Sedrakyan, and I. S. Burmistrov for discussions 
and I. E. Bulyzhenkov for valuable comments. The work 
was supported by RFBR, Programs of the RAS, Russian 
Ministry for Education and Science, and the Program 
"Leading Scientific Schools" . 



1 E. M. Lifshits and L. P. Pitaevskii, Statistical physics, part 
II: Theory of the condensed state. L. D. Landau course of 
theoretical physics, vol. IX (Pergamon Press, Oxford, New 
York, 1986). 

2 B. L. Altsfmler and A. G. Aronov, Electron-Electron In- 
teraction in Disordered Systems, edited by A. L. Efros and 
M. Pollak (North-Holland, Amsterdam, 1985). 

3 I. V. Gornyi and A. D. Mirlin, Phys. Rev. Lett. 90, 076801 
(2003). 

4 G. M. Minkov, A. V. Germanenko, O. E. Rut, 

A. A. Sherstobitov, V. A. Larionova, A. K. Bakarov, and 

B. N. Zvonkov, Phys. Rev. B 74, 045314 (2006). 

5 L. Li, Y. Y. Proskuryakov, A. K. Savchenko, E. H. Linfield, 
and D. A. Ritchie, Phys. Rev. Lett. 90, 076802 (2003) 

6 V. T. Renard, O. A. Tkacenko, Z. D. Kvon, E. B. Olshanet- 
sky, A. I. Toporov, J. C. Portal, Phys. Rev. B 72, 075313 
(2005). 

7 T. A. Sedrakyan and M. E. Raikh, Phys. Rev. Lett. 100, 
106806 (2008). 

8 V. G. Mokerov, B. K. Medvedev, V. M. Pudalov, 
D. A. Rinber, S. G. Semenchinsky, Yu. V. Slepnev, JETP 
Lett. 47, 71 (1988). [Pis'ma v ZhETF 47, 59 (1988)]. 

9 Throughout this paper the following notations are used: 
Ef = pp/2m* ,uj c — eB/m*,r — aora* /ne 2 , where m* is 
a band mass of the electron (0.067m e for GaAs and 0.21m e 
for Si) gd is the Drude conductivity. The latter is found by 
extrapolating the conductivity in ballistic regime to zero 
temperature. 

For short range scatterers case, r should be of the same or- 



der as the all-angle scattering time r q . From temperature 
dependence of the Shubnikov-de Haas oscillations ampli- 
tude we found that for our samples r q approximately equals 
r, indicating short-range scattering. 

A. Dmitriev, M. Dyakonov, and R. Jullien, Phys. Rev. 
Lett. 89, 266804 (2002); V. V. Chevianov, A. P. Dmitriev, 
and V. Yu. Kachorovskii, Phys. Rev. B 68, 201304(R) 
(2003). 

A. D. Mirlin, J. Wilke, F. Evers, D. G. Polyakov 
and P. Wolfle, Phys. Rev. Lett., 83, 2801, (1999); 
D. G. Polyakov, F. Evers, A. D. Mirlin and P. Wolfle, 
Phys. Rev. B 64, 205306 (2001). 

V. Renard, Z. D. Kvon, G. M. Gusev, J. C. Portal, Phys. 
Rev. B 70, 033303, (2004). 

G. Zala, B. N. Narozhny and I. L. Aleiner, Phys. Rev. B 
64, 214204 (2001). 

Because of additional valley degeneracy in silicon A = 1 + 
15F CT /(1 + F CT ) (T .A Sedrakyan, M. E. Raikh, private 
communication) . 

Fq value was determined from g-factor: V. M. Pudalov, 
M. E. Gershenson, H. Kojima, N. Butch, E. M. Dizhur, 
G. Brunthaler, A. Prinz and G. Bauer, Phys. Rev. 
Lett. 88, 196404 (2002). N. N. Klimov, D. A. Knyazev, 
O. E. Omel'yanovskii, V. M. Pudalov, H. Kojima, 
M. E. Gershenson, Phys. Rev. B 78, 195308 (2008). 
A. E. Voiskovskii, V. M. Pudalov, JETP Lett. 62, 947 
(1995). [Pis'ma v ZhETF 62, 929 (1995)]. 



